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Abstract 

The paper deals with the expected maxima of continuous Gaussian pro¬ 
cesses X = {Xt)t>o that are Holder continuous in L 2 -norm and/or satisfy the 
opposite inequality for the L 2 -norms of their increments. Examples of such 
processes include the fractional Brownian motion and some of its “relatives” 
(of which several examples are given in the paper). We establish upper and 
lower bounds for E maxo<f<i Xt and investigate the rate of convergence to that 
quantity of its discrete approximation E maxo<i<n Some further prop¬ 

erties of these two maxima are established in the special case of the fractional 
Brownian motion. 
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1 Introduction 

The problem this paper mostly deals with is hnding upper and lower bounds for the 
expected maximum 

E max Xt (1) 

0<t<r 

of a zero-mean continuous Gaussian process X = {Xt)t>o over a hnite interval [0, r]. 
In what follows, we will always assume that r = 1 which does not restrict generality 
in the case of deterministic r (note that the case of random r will require application 
of a more advanced technique, see e.g. [I9] and Section 1.10 in [H]). Computing 
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the value of expectation ([T]) is an important question arising in a number of applied 
problems, such as hnding the likely magnitude of the strongest earthquake to occur 
this century in a given region or the speed of the strongest wind gust a tall building 
has to withstand during its lifetime etc. 

We will mostly deal with processes X satisfying one or both inequalities in the 
following condition: for some positive constants Ci,Hi, i = 1,2, one has 

Ci\t-<\\Xt-Xsh <C 2 \t- for all t,s>0, (2) 

where ||f ^||2 = vE^ denotes the L 2 -norm of the random variable Y. 

In the case when Hi = H 2 in (2), the process X is a quasihelix in the respective 
L 2 space, in the terminology introduced in [13] (see also [Il|). So one can say that 
our paper mostly deals with the expected maxima of generalized quasihelices, with 
substantial attention paid to the special case of a helix hrst introduced in [15], which 
is the famous fractional Brownian motion (fBm) process with Hurst 

parameter H G (0,1]. This is a continuous zero-mean Gaussian process with = 0 
and covariance function 

EHf -\t- s\^^), t,s>0, (3) 

in which case \\B^ — B ^\\2 = \t — s\^. The Hurst parameter characterizes the nature 
of dependency of the increments of the fBm. For H G (0, |) and H G (|,1], the 
increments of B^ are respectively negatively and positively correlated, whereas in 
the special case H = ^ the process is the standard Brownian motion which has 
independent increments. When H = 1, the trajectories of the fBm process are a.s. 
rays with a random slope: B^ = ^t, t > 0, where ^ is a standard normal random 
variable. (It is instructive to note that, at the left end of the spectrum of H values, 
as hf —t 0, the limiting in the sense of convergence of hnite-dimensional distributions 
for B^ process will be just a Gaussian white noise plus a common random variable, 
see Section m) Other examples of processes satisfying ([2]) are given in Section (5] 

It is hardly possible to obtain a closed-form analytic expression for the expecta¬ 
tion dU), even in the classical case of the fBm. For the standard Brownian motion 
the exact value of the expected maximum is while for all other H 

within (0,1) no closed-form expressions for the expectation are known. 

One can try to evaluate expectation ([T]) numerically, using Monte Garlo simula¬ 
tions and the “discrete approximation” 

E max W ~ E max Xj/„, (4) 

0<i<l 0<i<n ' 

choosing a large enough n. One difficult question that arises in doing so is what 
value of n should be chosen to get a sufficiently good approximation. Fig. [T] presents 
the results of our attempts to estimate the value of ([T]) using approximation (|1|) for 
fBm’s with different Hurst parameter values. The twelve curves show the results for 
12 different values n = 2^,2^,..., 2^® (the values n = 2^ were chosen since paths of 
the fBm’s were generated using the Davies-Harte [1] algorithm which employs the 
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H 

Figure 1: Monte Carlo estimates for the values of Emaxi<j<„ n = 2^,..., 2^®. 
The fcth lowest curve shows results for n = k = 1,..., 12. 


fast discrete Fourier transform, and the latter is convenient to compute when the 
numbers of points is a power of two). For each combination of the values of H and 
n, we simulated 5 x 10® paths and the same number of antithetic paths. In all cases, 
the 99.9% conhdence intervals were of lengths less than 0.02. 

As seen from Fig. [H approximation (|1]) does not seem to work well for the values 
of H close to zero, as one can hardly notice in the respective part of the plot any 
convergence of the estimates to a particular value as n grows. 

This leads to the natural problem of bounding the difference 


A„(X) ;= E max — E max Xnn > 0 

0<t<l 0<i<n ' 


(5) 


for a Gaussian process X satisfying ([2]) and, in particular, for the fBm. For the 
latter, in view of the observed slow convergence, we are particularly interested in 
the case of small H values. 

The main results of the paper show that, for a process X satisfying ([2]), 


KiC, 

ViTi 


< E max Xt < 

0<t<l 


K2C2 


A„(X) < 


A'3C'2\/lnn 

nH2 


( 6 ) 


with some absolute constants iCj G (0, 00 ), z = 1, 2, 3. For a family of processes X^^, 
Hi G (0,1), satisfying the left inequality in (|2]) with Ci = Ci^ > Cq > 0 (or even 
Ci^ S> '/Hi), the hrst relation in ([H]) shows that Emaxo<t<i X/^ —)■ 00 as iJi —)■ 0, 
while the last one suggests a poor convergence rate for the discrete approximation. 

Furthermore, we provide some lower bounds for the discrete approximation rate 
in Section[3l In addition, using these bounds, we study the properties of the expected 
maximum of the fBm and its discrete approximation. We show that the latter is 
continuous as a function of H and hnd its limit as H —)• 0. 

Computing bounds for the extrema of Gaussian processes is a large research area, 
see e.g. monographs [301 ED and references therein. A powerful general method for 
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obtaining bounds for the expected maximum of a stochastic process is based on 
a generic chaining technique, for which [2^ ISP] are exhaustive sources. However, 
for the particular problem we are dealing with here, the general method seems to 
be unnecessary cumbersome (moreover, it includes a non-trivial step of choosing 
partitions of the time interval [30l Section 2.3]). The simple argument we provide 
below already allows one to obtain lower and upper bounds of the same order of 
magnitude 1 / \/H. 

There are quite a few results in the literature concerning a related problem on 
evaluating the so-called Pickands’ constant Ph that plays a fundamental role in the 
theory of extrema of Gaussian processes. The constant is dehned as 

Ph := lim -E max exp(V2Bf^ - 

r-foo r 0<t<T ^ ' 

Recent results on theoretical and numerical estimates for Ph include 0151191121]. In 
particular, m contains an auxiliary result giving an upper bound for the expected 
maximum of the fBm. 

With regard to the problem on bounding A„(X), note that there are many 
results about approximating the paths of continuous time processes by those of 
discrete time ones, that can yield bounds for their functionals, including maxima 
and minima. However, the path-wise approximation may not be necessary in the 
problem we consider, as we are only interested in the expected maximum. Regarding 
bounds for the expected maximum, we can only mention the preprint m, where 
the maximum of the geometric fBm was considered in the context of option pricing 
in the fractional Black-Scholes model. The bounds obtained in [TT] for the discrete 
approximation are of the order 0{n~^ yhrn), but the constant in the bound depends 
on H. 

The paper is organized as follows. In Section [2] we prove lower and upper bounds 
for the expected maximum. In Section [3] we establish an upper bound for A„(X) 
and also consider the problem of bounding the difference from below. In Section S] 
we establish properties of the expected maximum of the fBm and its approxi¬ 
mation as functions of H and n, in particular, their continuity and limits as if —)■ 0. 
Section [5] contains further examples of Gaussian processes satisfying ([2|), so that the 
results of our Theorems [TH5] are applicable to them as well. The Appendix contains 
some known inequalities for Gaussian processes that are used in the paper. 


2 Bounds for the expected maximum 

In this section we use a combination of Sudakov’s inequality and Fernique-Talagrand’s 
generic chaining bound to give short derivations of upper and lower bounds for the 
expected maximum of a centered Gaussian process satisfying ([2|). In particular, the 
bounds imply that, for a family of processes H G (0,1), such that satis- 
hes ([2|) with Hi = H 2 = H, the constants G* being independent of H, the expected 
maximum of on [0,1] is of the order of magnitude of 1/\/H as if —)■ 0, 
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Theorem 1. (i) If there exist C > 0 and H G (0,1) such that HX^—X 5 II 2 > C\t—s\^ 
for any t,s E [ 0 , 1 ] then 


E max Xt > - > —=. 

\/AH Tie In 2 5\/ll 


(ii) If there exist C > 0 and H G (0,1) such that \\Xt — Xs ||2 < C\t — for 
any t,s E [ 0 , 1 ] then 


E max Xt < LC 

0<t<l 


271 


H\n'^2 


erfc 



IQ2C 


(7) 


where L < 3.75 is the constant from the generic chaining hound (Proposition 0] in 
the Appendix) and erfc is the complementary error function. 


Remark 1. Note that the basic assnmption of continnity of X is not needed in 
part (i). In the general case, the bonnd will hold for EsnpQ<^<^Xt which can be 
defined as, say, in (0.1) in [22] • On the other hand, that assnmption is snperflnous 
in part (ii), where the continnity of X will follow from the assumed bound for 
||Xt — Xs|| 2 , the Kolmogorov-Chentsov theorem and the usual scaling argument for 
Gaussian random variables. 


Remark 2 . The lower bound from Theorem [T] is noticeably more accurate than the 
upper one. Fig. |2]presents the graph of the lower bound and the value of the expected 
maximum for fBm computed numerically by the approximation at n = 2^® points. 
The upper bound is far above these plots. A better upper bound for H < 0.5 was 
obtained in 12 a Lemma 6 ] using Borel’s inequality, but simulations show that it is 
still considerably greater than the true value m provides an additional term under 
the square root, but the proof there is more technical than the simple argument given 
below that nevertheless gives the same order of magnitude 1/\/R). For if > 1/2, a 
much better bound than ([7j) immediately follows from Sudakov-Fernique’s inequality 
(Proposition 131): 

E max Xt < GE max = C\f^. 

0<f<l o<t<i 


Proof. Without loss of generality one can assume that C = 1 both in (i) and (ii), 
for one can always switch to considering the process Xt/C. 

(i) Let n be the integer part of the value x := that maximizes the function 

x“'^-^log 2 X. Then, in view of Sudakov’s inequality (Propositional), 


E max Xt > E max Xj/„ > 

0<t<l 0<i<n ' 


log2(w+ 1) 
n‘^^271 


> 


jogRT 

X^^271 


1 

2\/ii7reln2 


0.2055... 

Vh 


( 8 ) 
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Figure 2: The lower bound l/(5\/^) for the expected maximum of the fBm’s for 
different values of if G (0,1) and the Monte Carlo estimates for that expectation 
using approximation (|1]) with n = 2^®. 


(ii) Consider the sets Tq := {1/2}, := ,j = 1,... ,2^"} for n > 1, so 

that \Tn\ = 2^”. Using the generic chaining bound (Proposition 0] in the Appendix), 
we have, setting a := 1 ln2 and changing variables u := that 


OO pQQ 

E max Xi < L V <L dx 

o<t<i L 

n=0 

rOG T poo 

= L exp{a(a; - He^^^)}dx <- 

Jo ® Ji 


= L 


Hln^2 V 2 Vifln=^2 


□ 


3 Bounds for the discrete approximation 


In this section we study the discrete approximation to the expected maximum of X. 
The main result is Theorem [2] giving an upper bound for the approximation error ([5]). 
A rather crude lower bound for the error is obtained in Theorem [31 

Theorem 2. Assume that \\Xt — Xs ||2 < C\t — for any t,s ^ [0,1] with some 
constants C > 0 and H G (0,1). Then, for any n > 2}!^, 


A„(X) < 


2C'\/h/n/ 4 0.0074 \ 

V (lnn)3/2/ 


Proof. As before, we will assume without loss of generality that C = 1. Suppose 
that n > 2}!^ is hxed and let 


efc := 2n ^^\/k\n.n, k >1. 
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From the continuity of X and monotone convergence theorem it follows that 


A„(X) = > E max X^/ nk+i — max Xj/, 

V 0<i<n'=+l 0<i<n'= 

k=l - - 

°° roo 

y / P( max Xj/„fc+i — max X^/^k >u)du 
^ ^ Jq V 0<i<n^+^ 0<i<n^ / 

oo 

=• y + Jk)- 


k=l 

oo 

E 

k=l 


’ Cfc 


k=l 


Clearly, Ik < e^. To estimate Jk, introduce the sets Tm := {j/m,j = 0,1,... ,m}, 
Uk{s) := [t G fc+i . t y'— s, ‘ 2 fik ^ ^ 2 n^ ^ and note that 


P( max X^/^k+i — max X^i^k >m) < / / P{Xt — Xs > 

V 0<i<n'=+l ' 0<i<n^ ' / 

teUkis) 


< 


0<i<n^ 

1 




E E 


\/^ , s u{2rJ)^ 

^ seT^k teUkis) ^ ^ 


exp 


u 


2(l/(2n'')) 


2H 


n ■ (n^)^ 1 2/0 k\ 2 H 

< - ’ —exp — u{2n) 

- 2H^u 2 ^ ^ 

where the second line follows from the well-known bound for the Mills’ ratio of the 
normal distribution: 

„2 


P(^ > a;) < 


a 




exp 


X 


X 


2a2 


, a; > 0, 


for ~ X(0, (T^), and the observation that Var (X* — X*) < {2rJ)~‘^^ for t G Uk{s), 
s G T^k. The inequality in the third line holds since |T„fc+i| = + 1 = n ■ + 1 

(note that not all the points t G T^k+i contribute to the sum). 

Now applying the bound a > 0, we get 


n ■ (n’^y / 1 2 / k\2H 

Jk < -—F=—exp( — eli2ryy^ 




n ■ n 


2^^+‘^y/^k In 


-exp(-22^+Mnn'’) = 
n 


n 


n 


fc(l-H-22«+l) 


23^+2^^ In yy 


k 


Setting z := n 


1 H 22 «^+i^ conclude that 


A„(X) < y efc +-^ 

23^+2y^f 


n z 

nn k 


(9) 


Bound the contribution to the sum on the right-hand side from the first summands 
in the terms: 

f;.» = 2Vi^f;^ = 2^j(„"). (10) 


k=l 


k=l 


n 


n 
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where g{x) := xLi_i/ 2 (l/a;), a; > 1, lAk^z) being the polylogarithm function. It is 
not hard to show that gi{x) := x{g{x) — l), a: > 1, is a decreasing function, gi{2) < 4, 
so that g{x) < 1 + 4/x for x >2. Therefore, the right-hand side of fllUp is less than 
2 n“'^\/lnn(l -|- 4n“^), n > 2}^^. 

Next note that since 2 ; = < 0.008 when n > 2}!^ and g 2 {x) := 

— ln(l — x)/x, X > 0, is an increasing function, one has YlT=i ~ < 

5 ( 2 ( 0 . 008 ) 2 ; < 1.005z. Hence the contribution to the sum from the second summands 
in the terms on the right-hand side of (jH]) does not exceed 


23H+2y/^lnn 23^+3 \/^(lnn) 3/2 

0.0074 2\/h^ 

^ (lnn)3/2 


2\/\\in 


where we used the inequality < 2 - 2 in 4 holds for 

> 2 since sup^^g^Q 2(1 — 2'^^)/H = —2 In4. 

The above bound, together with inequality (jO]) and our bound for the right-hand 
side of flTU]) , completes the proof of Theorem [21 □ 

Remark 3. In the special case when X = is an fBm, it is possible to obtain 
a similar bound in a simpler way, but the result will contain a constant depending 
on H: for any H G (0,1) and n > 1, 


An{B^) < n ^\/a{H) + Inn, 


where a{H) is a function of H such that a{H) —)■ cxd as i7 —)■ 0 (see flTT|) below). 

Let B^’’^ be a process whose trajectories are random polygons with nodes at the 
points {i/n, i = 0,1,... ,n. Then 

AniB^) < E max(Rf - Rf’”) = EF, 


where Y := maxi<i<„ Z- 


H,n 


Zf’" : = 


max (B 

{i—l)/n<t<i/n 


H 


B 


H,n^ 


z = l,2. 


being identically distributed random variables such that 


r^H,Tl 


d 

= n 


max (Bf 

0<t<l 



due to the self-similarity of the fBm. 

Therefore, for any A > 0, by Jensen’s inequality one has 

(EF)^ < -- In Ee^^^ = -- In E max ^ T ln(nEe^*'^i ’ ^^). 

A A l<i<n A ^ ^ 











Now taking A := we obtain the desired bonnd for An{B^) with 


a{H) := In Eexp 
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Again using Jensen’s inequality, we see that 
aiH) > E max (5f - )) > E (^max Bf - max{0, 

>E(maxSf) - E(^max{0, Sf}) >(EmaxBf^ - E(^max{0, , 


which tends to inhnity as —)■ 0 in view of Theorem [T](i). 

Remark 4. Note also that, in the even more special case of the standard Brownian 
motion (H = 1/2), the exact asymptotics of are well-known and contain 

no logarithmic factor: 

= (a -f- o(l))n“^'^^, n —)■ oo, 

where a = —C(l/2)/\/^ = 0.5826... and ({■) is the Riemann zeta function 
(see |Sl ES]). One may expect that the logarithmic factor in the general bound 
from Theorem [1] is also superfluous, but that the power factor n~^ gives the correct 
order of magnitude for An{B^). 

Unfortunately, in the general case we are not aware of any good lower bound for 
the discrete approximation rate. In the case of our particular interest, when if —)■ 0, 
our next theorem shows that the number of points n = n{H) should grow as an 
exponential of 1/if to provide a reasonable approximation. Namely, we will show 
that n(if) should be at least of the order of for some constant c > 1 or, which 
is the same, that {n{H))^ should be bounded away from below from one as if —)■ 0. 

Theorem 3. Let H G (0,1), be a family of zero-mean Gaussian 

processes such that, for some constants Ci, C 2 > 0 and all H G (0,1), one has 

Ci\t — s\^ < ||A/^ — < C' 2 |t — s|^ for all t,s & [0,1]. 


Suppose n{H) >2 is such that 

limsup A„(j:^)(A^) < 00 . (12) 

Then liminf//_j.o(n(if))'^ > 1. 

Proof. Suppose the contrary: relation flT^ holds, but there is a sequence if^ —t 0 
such that ffmlnn(fim) —t 0. 

Set Um ■= n{Hm) and choose km such that < log 2 n(ifm) < 2^"*. Note 

that, in view of the above assumption, one has 

2 *" = (13) 
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Now introduce the sets Tq := {1/2}, Tj := i = 1,..., 2^'’} for 1 < j < km, 

and := Tj,^ U {i/n^ : i = 0,1,...,n^}. Note that \Tj\ < 2^" for all j < fc^+i- 
Hence we can apply the generic chaining bound (Proposition 0]) with T = to 

get 


E max ^ + L 2 ^krr.+i)/ 2 (j 

0<i<nm ^ 

j=o 

r^m 


< LC 2 


/Hx + 


in view of fll3p . Setting a := |ln2 and making the change of variables u := e“^, we 
see that the integral on the right-hand side of the above formula equals 


exp[—aHme^°‘^)dx = 




2^rn,l2 

'du < -= 


again in view of fll3p . Thus, Emaxo<i<„^X/;(}^ = o{H, 
From here and Theorem [T](i), one has 


- 1 / 2 . 


Cl m-l 


An = E max — E max ^ _|_ o{Hn}C\ qq 

’ 0<t<l * 0<i<n™ - 5 ^ V m y 

as m ^ 00 , which contradicts (IT^ . Theorem Elis proved. 


□ 


4 Further properties of the expected maximum of 
fBm 

Consider the functions 

/(F, n) - E max B.%, //€(0,1). 

0 <i<l 0<2<n ' 

It follows from Sudakov-Fernique’s inequality (Proposition [1]) that both f{H) and 
f{H,n) are non-increasing in H. Theorem [T](i) implies that f{H) —)■ 00 as if —)■ 0. 

The main goals of this section are to show that f{H) and f{H, n) are continuous 
in H and hnd the limit of f{H, n) as ii —)■ 0. 

We will start with an auxiliary lemma about the limit of hnite dimensional 
distributions of , which is of independent interest. Let ^ be the standard 

Gaussian white noise, so that are i.i.d. standard normal random variables. Dehne 
the process Z = {Zt)t>o by setting Zt := (^ - ^ 0 )/\/2, t > 0. 

Lemma. As H ^ 0, Z in the sense of eonvergence of finite-dimensional 

distributions. 
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Proof. Since all the finite-dimensional distributions of and Z are zero-mean 
Gaussian, it suffices to establish convergence of the covariance function of to 
the one of Z. The latter is clearly given by 

f 0 if f = 0 or s = 0, 

^ZgZt = < I if f 7^ s and t,s > 0, 

1 if f = s > 0. 

That the covariance function ([3|) of B^ converges to the same expression as if —)■ 0 
is obvious. □ 


Theorem 4. (i) For any fixed n > 1, 


lim f(H, n) 
H^o ^ ^ 



where fi are i.i.d. standard normal random variables, := max{0 ,x}. 

(ii) For any fixed n > 1, the function f{F[,n) is continuous in H on [0,1). The 
function f{H) is continuous on (0,1). 

(hi) Let Hk G (0,1) and nk> 1 be two sequences, Hk 0 as k ^ oo. Then the 
sequence fk := f{Hk,nk) is bounded if and only if Uk is bounded. 

Note that part (iii) could easily be proved if we knew that the function f{n, if) 
is monotone in n. However, we do not have a simple proof of this fact, and therefore 
will present a proof based on the bounds obtained in the previous sections. 


Proof. Assertion (i) follows from the lemma, with the positive part of the maximum 
appearing because Zq = 0. 

(ii) Suppose 0 < ifi < ii2 < 1. As we noted above, f{Hi,n) > f{H 2 ,n) for any 
n > 1, and applying Chatterjee’s inequality (Proposition H]) we get 


0 < f{H,, n) - f{H 2 , n) < V«„(iii, fi2) Inn, 


where, for ifi > 0, 

aJHi,H 2 ) = max ((j/n)^-^^ — < max 

^<j<n 0<x<l 

H2 fjr, — ff 1 

= (iil/ii2)«2-«l - (ifi/if2)«2-«l < ^ , (14) 

eh I 

while, for ifi = 0, 

Q!„(0, H 2 ) = max (1 — (j/nY^^) = 1 — (15) 

l<j<n 

Now the continuity of f{H,n) at any H G (0,1) and at ii = 0 follows respectively 
from fflTD and flT^ . 

The continuity of /(ii) at any H G (0,1) is obvious from the continuity of 
f{H,n) and the upper bound for An{B^) = f{H,n) — f{H) from Theorem [2j 
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(iii) If rzfc is bounded, the boundedness of fk follows from (i) and (ii). In the 
case of unbounded sequence nk, it will suffice to prove that fk^oo provided that 
Hk ^ 0 and —)• cx) as fc —)■ cxo. 

It follows from Theorem [2] that, for n > 2^/^, one has 


/./rr X 6\/lnn 2-0.0074 6\/ln?7, 

f{H) - f{H, n) < -jT- + •= 0-0107. 


nH 2^/^ 


n 


Hence we obtain from Theorem [U^i) that 




evi 


nn 


5y/H 


n 


H 


-Cl, n>2 ^/^. 


On the other hand, in view of Sudakov’s inequality (Proposition |3]), for any n > 1 
one has (cf. ([8j)) 

f{H,n) > C2 := (27rln2)“^/^ 




Therefore, 


f{H, n) > max<| 

> 


6\/lnn C2\/lnn 


1 5^/H n 


C2 


(6 + C2)h\fH 


H ’ 

-Cl, n> 2^/-^, 


Cl 


and 


f{H,n) > ^Vlnn, n < 2^''^. 


Now the desired claim immediately follows from flT6|l and ffT7|l . 


(16) 

(17) 

□ 


5 Examples 

As we mentioned in the Introduction, a fundamental example of a zero-mean Gaus¬ 
sian process satisfying ([2]) is the fBm (in that case, ([2]) holds with Hi = H 2 = H, 
Cl = C 2 = 1). In this section we give some further examples of such processes. 

1. Sub-fractional Brownian motion. This is a centered Gaussian process = 
{Cl^)t>o with parameter H G (0,1), such that its covariance function equals 

EGfGf = + ^(|t + s|2^ + |t-s|2^), t,s > 0. 

This process was introduced in [T] in connection with the occupation time fluctu¬ 
ations of branching particle systems. In the case H = 1/2, it coincides with the 
standard Brownian motion: _ For H 7^ 1/2, is, in a sense, a pro¬ 

cess intermediate between the standard Brownian motion anij the fBm B^. 
Sub-fractional Brownian motion satishes the inequalities (see [1] for details) 

\t — s\^ < \\C^ — Cf \\2 < \/2 — 2^~^\t — s|^, t,s >0. 
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2. Bi-fractional Brownian motion. Let where H G (0,1), 

K G (0,1] are parameters, be a zero-mean Gaussian process with covariance function 

-\t- t,s>0. 

This process can be considered as an extension of the fBm, the latter being a special 
case when K = 1 (see uniEiE]). The process B^’^ satishes the following version 
of (ED: 

2-ir/2|^ _ < 2^i-K)/2^t _ ,\hk ^ ^ > q_ 

3. Fredholm and Volterra Ganssian processes. A large class of Gaussian 
processes can be obtained via a Fredholm representation: 

Xt= [ IC{t,s)dBl/^, tG[0,l], 

Jo 

with some kernel K, G -L^([0,1]^) and a standard Brownian motion 5^/^. In the case 
when /C(t, s) = 0 for all s > t, this is called a Volterra representation. 

It was proved in [SB] that a zero-mean Gaussian process on [0,1], which is sep¬ 
arable in the sense that the Hilbert space L‘^{Q,a{X),P) is separable (note that 
all continuous Gaussian processes are separable), admits a Fredholm representa¬ 
tion if and only if its covariation function R{t, s) := EXtXg satishes the condition 
R{t,t)dt < oo. Then it follows from Ito’s isometry that a sufficient condition 
for ([^ to hold is that 

Cl\t — < f {]C(t,u) — X{s,u))'^du < C^lt — for all t,sG[0 ,l]. 

Jo 

4. Wiener integrals with respect to fBms. For a non-random measurable 
function / : [0,1] —>■ R and a hxed H G (0,1), consider the process 

Xr.= f f{s)dB^, tG[0,l]. 

Jo 

If f{t) is Holder continuous with exponent a > 1 — H then the integral can be dehned 
as a pathwise Stieltjes integral. The thus dehned integral will coincide with the one 
given by the following dehnition employing fractional integrals which is applicable 
in the general case. 

Recall that the right-sided Riemann-Liouville integral {X^_f){t) of order a G 
[0,1] of the function / is dehned as follows: (I°_/)(t) := /(t), t G [0,1], while for 
a G (0,1] one puts 

(X"_/)(t) := f f{s){s - t)“-^ds, t G [0,1]. 

i («) Jt 
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For a G (—1,0), define (X“_/)(t) as the fractional derivative of / of order |a|: 

(ir_/)(() := 

Let be the linear operator that maps / : [0,1] —?■ R to the function 

(K^m) ( 6 [O.l], 

where (• denotes the function 1 1—)■ and Ch is the normalizing constant 

_/ 2Hr{3/2-H) y/" 

^ Vr(2-2R)r(R + i/2)y ’ 

The fBm admits the representation (see e.g. Section 1.8 in [T2] and ini) 

Bf = [\K«lpj,){u)dB'J^, «e|0,l], 

Jo 

where I[o,t) is the indicator function of [0,t). Now, for / such that f G X^([0,1]), 
the integral of / with respect to can be dehned as 

X, := [‘ f{u)dB« = [\K«lf I|„,„])(«)<iBy^ ( 6 |0,1], 

Jo Jo 

Next we will give simple sufficient conditions that ensure that relations (EJ with 
Hi = H 2 = H hold for the integral process X = (Xt)tg[o,i]- 

For H > 1/2, using Ito’s isometry, it is easy to see that a sufficient condition 
for the left inequality in ([2]) to hold is that either f{t) > c > 0 for all t G [0,1] or 
/(t) < c < 0 for all t G [0,1]. For the right inequality in (EJ it suffices that \f{t)\ is 
bounded on [0,1]. 

For H < 1/2, assume for simplicity that / is absolutely continuous on [0,1] and 
set h ■= H — 1/2. Then the fractional integral in the dehnition of exists and 

(X/_[(-)Vl[o,*)])(«) = f{t){t - u)%^t){u) - j f{v))'\io^t){v){v - ufdv 

J U 

= - u)\^t}{u) - h f - u^dv 

J U 

n'*/'(n)I[o,t)(n)(n - u)’^dv, 

which can be obtained by integrating by parts (see e.g. Section 2 in [23]). Then a 
sufficient condition for the left inequality in (EJ to hold is that, for some c > 0, 

f{t) > c, f{t) — h~^tf'{t) > c for all t G [0,1]. 
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Indeed, in this case, using the preceding formula one obtains that 

■ )Vl[o,p])(«) > c(Xt[( • )^I[o,)])(«) > 0 


and so, by Ito’s isometry, 

E|X, - - Bff. 

Jo 

Similarly, for the right inequality in ([2]) to hold it suffices that, for some c < oo, 
\fit) \ < c, \f{t) - h~hf{t)\ < c for all t e [0,1]. 


Then we have 

I(IM( ■ )'71[0,<)l)(«)l < 2c(lM(. ,])(«) 

and, again by Ito’s isometry. 


E\Xt - < 4c2E(5f - B^f. 

6 Appendix: Inequalities for Gaussian processes 

Here we collected a few classical bounds for the expectations of the maxima of 
Gaussian processes that we needed for deriving our results. 

Proposition 1 (Sudakov-Fernique’s inequality, [51 EZl El]). Let X = {Xt)o<t<i 
and Y = (F))o<t<i be two continuous Gaussian processes such that EXt = ET) and 
- Yt \\2 < IIX^ - Xt \\2 for all s,t ^ [0,1]. Then 

E max < E max Xt. 
o<t<i o<t<i 

Remark 5. Sudakov [27] and Fernique [^ proved the inequality from the claim 
of Proposition [T] under the assumption that the processes X and Y are zero-mean. 
Vitale [5T] extended the assumptions requiring only that EXt = EYt. 

Proposition 2 (Chatterjee’s inequality, [2])- Let (Xi,..., X„) and (Yi,..., Yn) he 
Gaussian vectors such that EXj = EF), i = l,...,n. Set Oij := ||Xj — XjWl and 
bij := IIV - YjWl, i,j = Then 

E max Xj — E max V < . / max laj,- — bif ■ Inn. 
l<i<n l<i<n - V l<i,j<n ' ^ ^' 

Proposition 3 (Sudakov’s inequality, [27]). Let (Xi,... ,X„) be a zero-mean Gaus¬ 
sian vector. Then 

E 


max Xi > y min ||Xi - Xj \\2 

l<i<n V 27r ij^j 
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Proposition 4 (Generic chaining bound, [28l |29l [30]). Let X = {Xt)o<t<i he a 
continuous zero-mean Gaussian process and Tq C Ti C T 2 C • • • C T C [0,1] a 
sequence of sets such that |To| = 1 and |T„| < 2^", n > 1. Then 



E max Xt < L max 2'^^‘^mm\\Xt- Xsh, L < 3.75. 


Remark 6. The literature on bounds for Gaussian processes usually does not pro¬ 
vide explicit bounds for the constant L. The bound L < 3.75 can be obtained from 
the proof in [29]. Note also that if the set T is hnite, then the sum in the bound 
will contain a hnite number of terms because in this case T„ = T for n large enough 
(excluding the trivial case when T„ “stop growing” at some point and one never has 


Tn = T). 
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